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We study the dynamics of a general quartic interaction Hamiltonian under the influence of dis¬ 
sipation and non-classical driving. We show that this scenario conld be realised with a cascaded 
superconducting cavity-qubit system in the strong dispersive regime in a setup similar to recent 
experiments. In the presence of dissipation, we find that an effective Hartree-type decoupling with 
a Fokker-Planck equation yields a good approximation. We find that the stationary state is approx¬ 
imately a squeezed vacuum, which is enhanced by the Q-factor of the cavity but conserved by the 
interaction. The qubit non-linearity, therefore, does not significantly influence the highly squeezed 
intracavity microwave field but, for a range of realistic parameters, enables characterisation of itin¬ 
erant squeezed fields. 


I. INTRODUCTION 

Open quantum systems methods are extremely valu¬ 
able for studying many systems of interest in quantum 
information and control, where the extent to which dissi¬ 
pation and decoherence can be reduced is limited by the 
need to pass signals in and out of the system [1]. These 
problems are particularly difhcult to analyse in the pres¬ 
ence of internal non-linearities and interactions [5] and 
possess no closed analytical solution in general, requiring 
a variety of approximate analytical and numerical tech¬ 
niques to proceed [5]. There is therefore great interest 
in gaining insight in situations in which the stationary 
state of the system is non-trivial yet can be analysed. In 
particular, the ability to produce, detect and characterise 
non-classical electromagnetic states is increasingly being 
explored [IHB]. 

In this paper we analyse the dynamics of a non-linear 
(quartic) open quantum oscillator which is driven by 
a non-classical field (squeezed vacuum [7]) and find its 
steady-state. The non-linearity is small compared to the 
dissipation and the drive is modelled by a cascade of 
another parametric oscillator and the non-linear oscil¬ 
lator. This problem is related to both the open and clas¬ 
sically driven Duffing oscillator [8] and the the problem 
of a two-level system interacting with a squeezed reser¬ 
voir [S], which are analytically solvable. Here, however, 
the model does not yield to these analytical techniques 
and we instead develop a combined Hartree and Fokker- 
Planck equation self-consistent treatment, which admits 
a class of Gaussian stationary states. 

This result is applicable to a variety of systems, but we 
focus on the case of a superconducting cavity-qubit sys¬ 
tem operating in the strong dispersive regime. We show 
that it is possible to generate a highly squeezed vacuum 
state in a superconducting resonator driven by the output 
of a Josephson parametric amplifier and, by comparison 
with simulations of the Lindblad master equation for the 
full, unsimplified system, show that our model describes 
the system well in this limit. In this parameter range, full 


state tomography of the cavity can be achieved using the 
qubit, giving the potential to use this system as a means 
of characterising a traveling squeezed field Hiiniiii]. 

Circuit quantum electrodynamics (cirQED) has pro¬ 
vided an excellent test-bed for fundamental quantum op¬ 
tics, owing to the largely dissipationless environment pro¬ 
vided by superconductivity |12j . The large non-linearity 
provided by the Josephson junction allows the produc¬ 
tion of high quality transmon qubits [nm] which can be 
strongly coupled to microwave resonators. One of the pri¬ 
mary benefits of cirQED is the ability to go beyond dis¬ 
persive quantum optics, where cavity-qubit detuning 
is much greater than their coupling g, and work with pa¬ 
rameters such that g^ jlS.q > k, the cavity width. In this 
strong dispersive regime, number-splitting of the cavity 
[3 HU [IS] allows full state reconstruction to be performed 
using high-fidelity qubit measurements [muH], provid¬ 
ing a valuable tool for quantum information processing 
Coherent driving of the cavity-qubit system can 
take advantage of dispersive cavity shifts to measure the 
qubit or map the qubit state to the cavity state [22j . 
In contrast, here we have a drive with zero mean coher¬ 
ence and finite squeezing, i.e. dominated by fluctuations, 
which requires a fundamentally different approach. 

Recent work has demonstrated that it is now possible 
to efficiently produce squeezing in a superconducting cir¬ 
cuit and study the interaction with a highly non-linear 
system which can be considered an effective qubit m- 
This experiment confirmed the prediction of Gardiner [5] 
that exposure to a broadband squeezed vacuum will mod¬ 
ify the T 2 coherence time of an atom depending on the 
axis of squeezing. Given these advances, it is a natu¬ 
ral to investigate what happens when when driving with 
a squeezed input in the opposite limit, that of a weak 
non-linearity, which arises when the qubit is far-detuned 
from the cavity resonance. The interaction between a 
squeezed state and an on-resonance qubit in a cavity has 
been studied by Milburn in a closed-system context. 
Recently, squeezing in cirQED, in a similar setup to ours, 
has been proposed as a means to improve quantum state 
measurement as compared to coherent driving [24) . 
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The development of high quality Josephson paramet¬ 
ric amplifiers (JPAs) [25] has been vital to these devel¬ 
opments, and is still an active area of research. There is 
therefore a need for a good characterisation procedure to 
compare these devices. Homodyne detection methods are 
more challenging to implement in superconducting cir¬ 
cuits than conventional optics, requiring additional JPAs 
to amplify the signal and introducing additional noise [Bj . 
It is also unclear whether any distortions in the observed 
state originate in the source or the measurement am¬ 
plifier. Wigner tomography using a cavity-qubit there¬ 
fore has the potential to produce higher-hdelity measure¬ 
ments of an incoming squeezed field, while also providing 
information about non-idealities, contained in higher mo¬ 
ments of the field. 

We describe how we construct a Gaussian mean-field 
model of a non-linear system driven with the squeezed 
output of a parametric amplifier in Sec. |ll| and we go 
on to describe how this system could be implemented 
in cirQED in Sec. |III| Finally, in Secs. |IV| and |V| we 
discuss solutions of our model and compare these results 
with numerical solutions of the full quantum system. 


II. THE QUARTIC OSCILLATOR MODEL 

We begin by making a key observation regarding the 
closed part of the system dynamics, which is described 
by the Hamiltonian 


H = Loa'a + aa'a, (1) 

where a is the bosonic annihilation operator and w and C 
describe the resonator frequency and interaction strength 
respectively. For the class of squeezed vacuum states, 
(a) = 0 and, therefore, a simple mean-held treatment of 
this system will yield only trivial dynamics. Instead we 
wish to approximate H by some self-consistent Hamil¬ 
tonian depending on second moments of the cavity op¬ 
erators. As the uncertainty associated with a is given 
by Aa = (aa) — (a)^ = (aa), these moments represent 
Gaussian huctuations around the zero mean. Over suf- 
hciently short timescales (or, importantly, in the open 
system case, when is small compared with dissipation) 
and with an initial state that is at least approximately 
Gaussian, we expect that these terms will dominate the 
dynamics of the system. We therefore apply a bosonic 
Hartree-type approximation |26] to the Hamiltonian of 
the oscillator and obtain the second order Hamiltonian 

Heff = oja^a + ({aa)a'a' + C(a^a^)oo, (2) 

where u) = a;-|-^(4(a'^a) -I- 1) and we have neglected addi¬ 
tional terms that contain no operators and therefore do 
not effect the equations of motion. Naively, this appears 
to be a dentuned parametric driving Hamiltonian which 
should produce squeezing, but we can write down the 


Heisenberg equations of motion 

— (o'a) = 2i(^{a^a^)aa — 2i(^{aa)a^o' (3) 

at 

— (aa) = —2iCjaa — 2iC^{aa)(2a)a 1), (4) 

at 

and therefore show that 

^(a^a)=0 (5) 

— (aa) = —2i (2C(a^a) -I- C + w) (aa). (6) 

The average number of photons in the cavity remains 
constant, and consequently (aa) undergoes purely phase 
evolution. We also see that if we select Cj = —C(2(a'^a) -I- 
1) then we can achieve a stationary state. Note that our 
approximation includes the assumption that we can write 

(a'^oo'^a) « (a'^oo'^a)ii; = 2{a)a)"^ + |(aa)p -I- (a'^a), (7) 

and the model will break down if this reduced form is 
too different from the exact expectation value. However, 
we find there is a significant region of parameter space 
where this is not the case, which we show in Fig. and 
in our chosen application in cirQFD the model is valid in 
a regime where significant intracavity squeezing can be 
achieved. 

We now wish to study the behaviour of this system 
when driven with squeezed vacuum. It is well known that 
a parametrically driven resonator cannot achieve squeez¬ 
ing of more than a factor of two pTIEB] so instead we 
pump Hgff with the output of a parametric amplifier, 
described by 

Hi = wiojoi -I- -(eioj^ — e^af), (8) 

where wi is the frequency of the resonator and e encodes 
the drive strength and phase. The subsystems are con¬ 
nected by a uni-directional dissipative channel that con¬ 
nects the output of the ampliher to the input of second 
cavity, and the combined system is coupled to a zero- 
temperature bath. A formalism has been developed to 
study such open systems . which we use to obtain 

the full Hamiltonian of our two cavity system 

H = Hi+ Heff - (a^a'^ - a|a). (9) 

This is coupled to the bath by the combined collapse 
operator C = ^/kTui + ^/ku, where ki,k are the de¬ 
cay constants for the two sub-systems. The combina¬ 
tion of anti-symmetric coupling term in the sub-system 
operators and symmetric collapse operator results in uni¬ 
directional coupling between the two cavities. By moving 
into a rotating frame defined by uJia\ai we simplify the 
system to 


H = -{eitti — e*aj) -I- Aa^c 


I ^ / 12 * 2\ 

+ 21 “ - e a ) - 


iy/1^ (aio''’ - a\a), (10) 
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where we have defined A = w — wi and e = —Ci{aa). 

The standard quantum optics method is now to solve 
the master equation p = —i[H,p\ + C[C]p numerically, 
where C[C] is the Lindblad superoperator associated with 
the collapse operator C. Instead, we cast the system in a 
Fokker-Planck equation in the complex P-representation 

m. 


m 


P{a) 


d Idd 

oui 2 oai oaj 


P{a), 


( 11 ) 


where Aij and Dij are the (constant) drift and diffusion 
matrices respectively. The internal steady state spectral 
matrix (Fourier transformed covariance matrix) can be 
expressed in terms of these matrices using the relation 

S{uj) = ^{A + iujiy^D{A'^-iujiy\ ( 12 ) 


The integrated matrix Sij is a 4 x 4 object containing 
all possible second moments of the system. From this 
we can find the width of the cavity state in an arbitrary 
direction. For example the uncertainty in the quadrature 
P = - a) is given by AP = -S 33 - S '44 + 534 + 

^43 + 1/2. For our system the drift matrix is 


A = 


Kl , 

Cl 

* Kl 

*=1 2 

— ^KiK 0 

V 0 -i/k/K 


0 0 \ 

0 0 

-f-iA e 

e* * A - f / 


(13) 


while D = diag(ei, 63 , e, e*). A naive analysis of this lin¬ 
earised system, by looking at the real part of the eigenval¬ 
ues of A, suggests that this system possesses two thresh¬ 
olds where the stability of the fixed point at (a) = 0 
changes. The first is the well-known threshold of the 
parametric amplifier at ei = ki/ 2, with a second at 
\/kP-A 2 = k/2. In practice, however, reaching this 
threshold would require conditions that cause the Gaus¬ 
sian approximation to break down, a point which we 
briefly expand on in Sec. m 

Clearly, the value of e is not a free parameter and is in 
fact determined by the other parameters of the system via 
the value of the (aa) correlation function. This value is 
given by the entry S 33 of the integrated spectral matrix. 
To enable us to calculate the squeezing in the cavity, we 
evaluate e self-consistently, substituting the value back 
into A and D until the value converges. 


III. REALISATION IN CIRCUIT-QED 

Our superconducting cavity-qubit system is described 
by the Jaynes-Cummings Hamiltonian 

iJ 2 = ^ 20^02 + + 3 ( 020 -■*■+a|(T“), (14) 

where UJ 2 is the cavity frequency, g is the cavity-qubit 
coupling, ujq is the qubit transition frequency and are 


the qubit raising and lowering operators. The source of 
squeezing is a Josephson parametric amplifier pumped 
at twice the cavity frequency, which is described well by 
Hi [33]. In the strong dispersive regime, the qubit is 
sufficiently far detuned from the cavity that, provided 
the number of photons in the cavity remains low, it is 
never significantly excited. This allows us to eliminate 
the qubit by first diagonalising it block-wise to some or¬ 
der in a small parameter gjAq, where Aq = uiq — uji is the 
qubit-cavity detuning. Details of this procedure can be 
found in Ref. [T]. Here we take terms up to 0{g'^/A^^), 
leaving 


H 2 = {uj2-0ala2+u}qY+X 


4^2 + 1 0-2-5(0^02)^02, 


(15) 

where x = 9 '^ l\- 4/A, and 5 = g'^/A^, and set 0-2 = 
— 1. An additional effect of the diagonalisation procedure 
is to change the interaction between the two sub-systems 
to 


Hi = - 


^4^^(oiaJ-ola2) (^1+^02^ , (16) 


while Hi is unaffected. The combined system collapse 
operator is also transformed to 


C = i/iyai + v^02 ■ (17) 

As H 2 is of the same form as H^ff, it can be treated 
with the same Hartree approximation and so the super¬ 
conducting system of interest is described be our model 
with o —>■ 02 , e —>■ £2 = — 215 ( 0202 ), K —>■ /«2 = 
^ 2(1 — g^j and A —> A 12 = a ;2 — wi — x + 25 ( 0 ^ 02 ). 

In addition to out analytic results, we can construct 
a master equation for the full cascaded system including 
all terms in J/ 2 - We use the Qutip library [311 obtain 
expectation values of observables and reconstruct Wigner 
functions for the separate cavities as a function of time. 


IV. ENHANCED INTRA-CAVITY SQUEEZING 

As we show in Fig. this set-up can achieve up to 7 
dB squeezing of the intra-cavity field in the presence of 
a qubit and for parameters such that dispersive shifts of 
the cavity would allow full state tomography of the cav¬ 
ity. To demonstrate this, we consider the case of a high 
Q second cavity, also considered by Collett and Gardiner 
E7], where ki ^ K 2 - Specifically we take ki/k 2 = 50 
where ^2 = 1 MHz. In this case the squeezing is effec¬ 
tively infinitely broad compared with the cavity that is 
being driven. In the absence of the qubit, 3 = 0 and 
system can be solved exactly to reproduce their results. 
The best squeezing for any e is found at A 12 = 0 and 
AP tends asymptotically to 0 as ei approaches thresh¬ 
old. Strikingly, when the qubit is introduced we see that 
an identical degree of squeezing can be achieved, but 
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FIG. 1. (Colour Online) Plot of the P-quadrature uncertainty 
in the second cavity obtained in our theoretical model as a 
function of the effective cavity-cavity detuning A 12 and the 
pump strength ei. The other system parameters are fixed at 
ki/ac 2 = 50, sr/ft 2 = 56, Aq/ k .2 = 600. The total squeezing in¬ 
creases with pump strength, reaching the same uncertainty as 
the no-qubit case for any value of ei . The optimum detuning 
shifts as the total squeezing in the cavity is increased. 


at a non-zero value of A 12 , which grows with ei. This 
squeezing-dependent shift is shown in Fig. and occurs 
in addition to the number-dependent shift. As ei —>■ 0 
the position of this minimum tends to A 12 = —be¬ 
having like the solution of the undriven, dissipationless 
model in Eqn. in the limit of weak driving. At the op¬ 
timum Ai 2 for each ei, the axis of squeezing is always 
aligned with the incoming field, rather than at an angle 
as occurs with detuned linear cavities. This further sug¬ 
gest that the effective frequency of the cavity has shifted. 
This shifting, combined with the tightening of the ’val¬ 
ley’ in which any squeezing is seen, prevents the apparent 
threshold at ^|ep — A^ = /v/2, seen above, from ever be¬ 
ing reached as, while e increases with greater squeezing, 
so does the optimum value of A at which this squeez¬ 
ing is seen. To reach this threshold would require a very 
large qubit non-linearity, in which case our Gaussian and 
dispersive approximations would break down. 

As the effect of the qubit in this approximations is 
merely to shift the cavity resonance, we are able to pro¬ 
duce and observe squeezing much greater than a factor 
of two that can be achieved for the internal field of a 
parametrically driven cavity. The ability to produce high 
quality squeezed electromagnetic states is a valuable re¬ 
source for applications in high precision measurements 
of weak signals such as gravitational waves |35| : devel¬ 
opment of higher signal-to-noise communication proto¬ 
cols [36]; and provides a source of entangled photons for 
quantum technology such as key distribution [37j . At 
optical wavelengths, squeezing of 12.7dB below vacuum 



FIG. 2. (Golour Online) Plot of the minimum quadrature 
uncertainty that can be obtained in our system with K-ijK 2 = 
50 and g/K 2 = 56 by choosing the optimum value of A 12 (i.e 
the base of the ‘valley’ in Fig. [^. Mean field values in the 
Hartree approximation, which are identical to the no-qubit 
case, represent ideal squeezed states and fall on the boundary 
of the shaded region, which is inaccessible. Yellow circles 
show simulation data for no-qubit case, with deviations from 
theory caused by Fock basis truncation. Red diamonds and 
green squares show simulation data with Aq/ k .2 = 1200 and 
^qli ^2 = 600 respectively. With the qubit detuned further 
from the cavity, higher order interaction terms are smaller, the 
qubit behaves more like a spectator and greater squeezing can 
be observed. With Aq/K 2 = 1200 and t-ijK 2 = 12 we find 7dB 
of squeezing can be achieved, much greater than the factor of 
two | 28[ I29| that can be achieved in the internal field of a 
degenerate parametric amplifier, here marked with a dashed 
line. 


noise can now be achieved in a beam [SS] , but squeezing of 
an intra-cavity field has not been directly measured. The 
production of these states has been the subject of much 
recent work, investigating methods such as modulation of 
the cavity decay rate |S] , fast switching of the cavity res¬ 
onance |39j , and using parametric resonance driving [40] . 
New methods of state reconstruction, such as by sideband 
spectroscopy of the qubit [nma, have also been devel¬ 
oped. The ‘distillation’ of squeezing that we see has been 
discussed for a linear cavity driven with squeezed vacuum 
but it is not obvious that it should survive the qubit 
non-linearity. We are not aware of any reports of signifi¬ 
cant intra-cavity squeezing in experiment, but the level of 
squeezing we see is similar to that in recent experiments 
for itinerant squeezed states in superconducting circuits 
PiSllll] and is achieved in a simper set-up than other 
theoretical discussions, requiring no time-dependent pa¬ 
rameters. 


V. NON-GAUSSIAN STATIONARY STATES 

We compare our theoretical squeezing values with the 
results of full numerical simulations of the system master 
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(a) (b) 



(c) (d) 




FIG. 3. (Colour Online) Wigner functions of the squeezed 
cavity state in the second cavity with parameters hxed at 
ki/ac 2 = 50, p/k 2 = 56, Aq/K 2 = 600 and using the optimum 
value of Ai 2 (a) In the no qubit case with ci/ki = 10 the 
state produced is purely Gaussian and significant squeezed, 
with artifacts caused by basis truncation, (b) If the full qubit 
interaction is included rippling can be seen in the Wigner 
function, which damages the squeezing, (c) If only terms 
up to fourth order are considered, almost identical rippling is 
seen, suggesting that the breakdown of our mean field approx¬ 
imation is responsible for such distortions, (d) If the pump 
is increased to ei /«;2 = 12 the distortions become larger and 
damage the squeezing further. 

equation, 

P = -i\Hi + H 2 - ^ (ai4 - a\a2),p] 

-I- + y^a2]p, (18) 

both with the qubit and without (by setting ujq,g = 0) 
in Fig. This system includes all orders of the qubit 
non-linearity and therefore allows us t test the validity 
of our model. As the first cavity is fast and accumulates 
very few photons, we consider only 10 basis states while 
using 50 basis states for second cavity. We see that the 
simulations for the no-qubit case agree exactly with the¬ 
ory up to Cl = 13 (52% of threshold), above which there 
is significant deviation due to the truncation of the Fock 
basis. We only run simulations including the qubit for 
pumps below this value and test two different qubit con¬ 
figurations, both satisfying the number splitting criterion 
but with one qubit twice as far detuned from the cavity. 
We sweep over A 12 , with all other parameters fixed, to 
find the maximum squeezing at each the pump strength. 

With g/Aq « 0.1 we see very strong deviation from 
the model, preventing uncertainties of less than 0.2 from 
being achieved and fluctuations greater than the vacuum 



FIG. 4. (Colour Online) Plot of percentage error between the 
factorised fourth-moment {a^aa^a)R and the exact moment 
{a^aa^a). Green squares correspond to = 600 and red 

diamonds correspond to A 5 /K 2 = 1200 , matching the curves 
in FiglU with other parameters set to (ji/K 2 = 56, ki/k 2 = 
50. For the closer qubit, the error increases at relatively low 
pump strengths and is more than 5% for the largest values 
of e. This corresponding to the destruction of squeezing and 
non-Gaussian steady-state Wigner functions we see for these 
pumps in full simulations. For the further detuned qubit, 
however, the error is very low up to ei/K 2 = 13, leading to 
much better agreement between the model and simulations 
when Aq/ k .2 = 1200 . 

for high pump strengths. However, the optimal values 
of Ai 2 from the model are reproduced. In contrast to 
both coherent driving of a quartic interaction and uni¬ 
tary evolution under the quartic interaction, the system 
does attain a steady state and by plotting its Wigner 
function, as shown in Fig. we can see that the inter¬ 
action with the qubit introduces significant distortions, 
which increases with pump strength. This can be at¬ 
tributed to the breakdown of the Hartree approximation 
as the non-Gaussian part of {a^aa^a) term grows. In Fig. 
1^ we plot the percentage difference between the fourth 
moments and their factorised form from the model, and 
see a corresponding growth in the this error as we would 
expect. 

Doubling A, reduces the size of these higher order 
terms and the range of ei over which there is agree¬ 
ment with the model greatly increases. We see equiv¬ 
alent drops in the percentage error in the factorised mo¬ 
ments and the size of the distortions in the steady-state 
Wigner function. For these parameters, it is feasible to 
produce a highly squeezed state with AP « 0.1 and, as 
g^!A K, 2.6 > K, Wigner tomography can be performed 
experimentally. Additionally we see that the number dis¬ 
tribution of these states is dominated by even Fock states. 
This signature of pure squeezed states has not yet been 
observed in experiment and is illustrated in Fig. 

In order to use this set-up to perform accurate charac¬ 
terisation of itinerant squeezed vacuum, it is necessary 
to consider more higher squeezed incoming fields. In 
this case, it may be more appropriate to consider the 
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N 

FIG. 5. (Colour Online) Probability P{N) of observing N 
photons in the second cavity in the steady state. Solid bars 
show simulated number distribution in the steady state with 
K 1 /AC 2 = 50, ei/K 2 = 10, g/K ,2 = 56, Aq/K 2 = 600 and 
^i 2 Ik 2 = 4.95. Empty outlines show a comparison with an 
ideal squeezed cavity state of the same average photon num¬ 
ber. Almost all additional probability is in the N — 0 state, 
which is not shown to make these differences clearer. The 
simulated distribution is in good agreement with the ideal 
even-odd behaviour, which has not been observed experimen¬ 
tally to date. 


limit Ki « K 2 , to reduce distortions to the internal field. 
A good understanding would also be required of how 
higher-order non-linearities in the source effect the re¬ 
constructed field. A detector of this type would be of 


great utility in experiments and we plan to develop a 
more general model of squeezed driving to allow us to 
consider more highly squeezed and imperfect inputs. 

VI. CONCLUSIONS 

In conclusion, we have developed an effective model 
for a driven-dissipative system undergoing a quartic in¬ 
teraction and shown it possesses approximate squeezed 
vacuum stationary states under parametric driving. We 
have shown how this model arises in strong dispersive 
circuit-QED and that, in this setup, it is possible to gen¬ 
erate a significant intracavity squeezing in the presence 
of a qubit. We predict greater squeezing than has been 
achieved before in such systems, in a range of experimen¬ 
tally accessible parameters where dispersive cavity shifts 
enable state reconstruction. These results have potential 
application in the characterisation of sources of itinerant 
squeezed fields in superconducting circuits. 
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